Quantum Critical Behavior of the Cluster Glass Phase 
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In disordered itinerant magnets with arbitrary symmetry of the order parameter, the conventional 
quantum critical point between the ordered phase and the paramagnetic Fermi-liquid (PMFL) is 
destroyed due to the formation of an intervening cluster glass (CG) phase. In this Letter we discuss 
the quantum critical behavior at the CG-PMFL transition for systems with continuous symmetry. 
We show that fluctuations due to quantum Griffiths anomalies induce a first-order transition from the 
PMFL at T=0, while at higher temperatures a conventional continuous transition is restored. This 
behavior is a generic consequence of enhanced non-Ohmic dissipation caused by a broad distribution 
of energy scales within any quantum Griffiths phase in itinerant systems. 



In magnetic systems with disorder, it is always pos- 
sible to find isolated, defect-free regions which act as 
finite-size copies of a clean, bulk sample. Fluctuations 
of these rare regions (or droplets) lead to singularities of 
the free energy in what is known as the Griffiths phase 
[lj]. While these singularities are extremely weak (es- 
sential) singularities in classical systems [2j, the situa- 
tion can change dramatically near (T=0) quantum phase 
transitions 

B S, S, S S B 11- Recently, there has been 
significant progress made in understanding the quantum 
Griffiths phase by focusing on the dynamics of the rare 
regions jlOj ] - By considering the rare regions as inde- 
pendent droplets described by a local order parameter, 
an elegant classification of quantum Griffiths phenomena 
was proposed based only on general symmetry principles 
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In the independent droplet picture, the system is as- 
sumed to be correlated over the length scale of the in- 
dividual droplets, each of which has dynamics in the 
imaginary time direction related to the energy cost of 
coherently flipping a large volume of spins. The lo- 
cal droplet degrees of freedom therefore map to clas- 
sical, one-dimensional ferromagnetic spin chains in the 
imaginary-time direction, with their dynamics corre- 
sponding to interactions in this single dimension. The 
critical behavior of the bulk system can then be un- 
derstood through the behavior of the independent spin 
chains which depends only on their dynamics and the 
spin symmetry For example, in the case of itin- 

erant Ising magnets, independent spin chains are effec- 
tively above their lower critical dimension and undergo a 
phase transition in the universality class of the Kosterlitz- 
Thouless transition Thus, if sufficiently large, these 
droplets freeze over, leading to their eventual ordering 
at low enough temperature. Therefore, at zero tem- 
perature, the system is always in its ordered phase and 
the quantum phase transition present in the clean mag- 
net becomes destroyed or rounded due to the presence 
of disorder 0). In contrast, independent droplets in 
itinerant Heisenberg magnets are exactly at their lower 
critical dimension corresponding to weaker dissipa- 
tion. Within the independent droplet picture, this be- 



havior leads to quantum Griffiths effects, and the quan- 
tum phase transition to the magnetically ordered phase 
is conjectured to be in the universality class of the infinite 
randomness fixed point (IRFP) Q. These general issues 
have attracted much interest in recent years due to the 
possible role of such quantum Griffiths phase phenomena 
as a driving force for the disorder-driven non-Fermi liq- 
uid behavior in correlated electron systems (for a recent 
review see Ref . IT3 ) . 

In realistic systems, however, the droplets are not in- 
dependent and interactions must be considered, provid- 
ing a qualitatively different mechanism for dissipation. 
Such interaction effects were recently [j| shown to rep- 
resent a singular perturbation within any quantum Grif- 
fiths phase, leading to enhanced non-Ohmic dissipation. 
This mechanism leads to the freezing of sufficiently large 
droplets and the generic emergence of a cluster glass 
(CG) phase intervening between the uniformly ordered 
phase and the paramagnetic Fermi liquid (PMFL). (For 
a schematic phase diagram of an itinerant magnet with 
and without interactions, we refer the reader to Ref. |tt) 

In this Letter, we examine the nature of the quantum 
phase transition to the cluster glass phase. To study the 
critical behavior at the quantum CG-PMFL transition, 
we focus on the continuous spin symmetry case for which 
the interaction effects produce a sharp phase transition 
at T = 0. By studying the transition both from the 
CG and from the PMFL, we show that a region of co- 
existence of these two phases exists, indicating that the 
T = quantum phase transition is first-order. At higher 
temperatures, a tricritical point is found above which a 
conventional second-order transition is restored. In con- 
trast, we show that a collection of identical droplets dis- 
plays a conventional (continuous) phase transition even 
down to the lowest temperatures. From this we conclude 
that fluctuations due to Griffiths phase anomalies are re- 
sponsible for driving the transition first-order. 

Self- consistency conditions. — Our approach here mir- 
rors that of previous works [8, 9]. In the presence of dis- 
order, defect-free regions can form with near uniform or- 
dering while the bulk system is in the non-ordered phase. 
These rare regions (droplets) can be expressed as a flue- 
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tuating order parameter field in the spirit of Ginzburg- 
Landau with effective action S = ^ ^L.i+X^j ■ For 
the case of an itinerant antiferromagnet [13( , the local ac- 
tion (for a single droplet) is 

= J~] <f>i(u n )(n + \u n \)<t>i(-w n ) (1) 



where </>i is the iV-component order parameter field for 
the i th droplet and j-j is its bare mass, selected from the 
distribution Viji). This local action has been studied 
previously and is well understood |8j. The novelty arises 
when we consider interactions between droplets, medi- 
ated by the weakly correlated metal in the bulk. These 
take the form of the RKKY interaction 



J j > 







dr(pi(r)<t)j(T), 



(2) 



with Jij random, of zero mean and having variance 
(Jfj) = J 2 - Very recent work has suggested that even 
infinitesimally weak RKKY interactions destabilize the 
Griffiths phase found in the non-interacting droplet the- 
ory, leading to the cluster glass phase [9|. These results 
demonstrated that the nontrivial physics of non-Ohmic 
dissipation persists even when the problem is solved at 
the saddle-point level, which is formally justified in the 
large- N limit of the model [l4[ . 

We proceed in the usual way by averaging over disor- 
der using the replica method and decoupling the resul- 
tant quartic term by introducing an auxiliary Hubbard- 
Stratonovich field [jj, [HI]. This new field describes the 
dynamics of spin fluctuations caused by the long-ranged 
RKKY interactions, and contributes an additional dissi- 
pative term to the action. Its saddle-point value is [9( 
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dn 



n + Xi + \u n \ - gx{u n ) 



(3) 



where g = J 2 ~ 2<i is the RKKY coupling. A, 

is a Hubbard-Stratonovich field decoupling the quartic 
interaction in the single-site action and satisfies the local 
self-consistency condition 14 1 



A, = — 



uT x 



1 



n + \i + \cj n \ - gxi^n) 



(4) 



Our original action now becomes a single-site problem 
with effective action 

S c g = ^2 (f>i(oJ n )(ri + Xi + \<jj n \ - gx(u n ))4>i(uJn)- (5) 



Physically, this corresponds to independent droplets in 
the presence of additional dissipation arising due to in- 
teractions with all other droplet degrees of freedom and 



represented by the u> dependence of the local spin sus- 
ceptibility x(w„). 

Quantum criticality for uniform droplets. — An inter- 
esting limit of this problem is revealed if all the droplets 
are assumed to have the same size V(ri) — S(ri — r) for 
all sites i. In this case, the coupled set of equations © 
and ([4^_ are identical to those describing a metallic spin 
glass [16( and, at large enough g, we expect a zero tem- 
perature transition from a paramagnet to a spin glass. 
The self-consistency condition ([3j) for x(w) reduces to a 
simple algebraic equation and is easily solved: 



(6) 



±V(r + \+ \uj, n 



2ff 



where A = A(f). Clearly, the paramagnetic solution be- 
comes unstable when (f + A) = \/2g and we define an 
energy scale 17| A = (f + A — \/2g) . The phase bound- 
ary can be determined by setting A = and solving the 
self-consistency condition |4]) for r(T). 

The phase for A > (f > f c ) is characterized by the 
frequency dependence of Xu(w n ). From the solution ©, 
we find three separate regimes, 



Xu(o) - xuM 




,w < A < ^/2g 
, A < uj < y/2g 
,w> y/2g 



(7) 



The linear, low-w behavior of % u is characteristic of a 
Fermi liquid, while, right at the transition, % u ~ ^/uj, 
characteristic of non-Fermi liquid behavior. 

The transition at T = can also be studied from the 
magnetically ordered side. In this case, the quantum 
critical point is identified when the mean-field stability 
criterion for the glass phase vanishes [HI i.e. Asg = 
1 — \/2~9x(0) = 0- Again, using the solution ([6]) for \, we 
find Asg ~ vA which vanishes at A = 0. 

Thus, examining the transition through the instability 
of either the paramagnetic Fermi liquid or the spin glass 
phase yields a quantum critical point at A = 0, consistent 
with a conventional second-order transition in this limit 
of the problem. 

In the following we consider the model with a distribu- 
tion of droplet sizes, and show that the two approaches 
do not coincide at the transition, revealing the singular 
effects of droplet fluctuations. 

Griffiths phase behavior of distributed droplets. — For 
dilute impurities, the defect-free regions assume a Pois- 
son distribution, which can be written in terms of the 
local coupling constant as 



2-rra(ri —r) j u 



r,<f< 0. 



The offset f is used to tune through the transition, a is 
related to the droplet density p by a = pu/2it6 and <5 < 
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is the coupling constant of the clean magnet. In this way, 
the distribution is smoothly connected to the uniform-r^ 
limit as a — ► oo. We can therefore reasonably apply 
the uniform solution ([5]) as a zeroth order approximation 
and iterate equations and (0| until self-consistency 
is achieved. By accounting for droplets of all sizes, we 
are also including their fluctuations which are associated 
with the Griffiths behavior. 

The instability of the Griffiths phase at T = is al- 
ready apparent at one iteration loop (IIL). Physically, 
this corresponds to the regime where very large droplets 
are very dilute, but the typical droplets (in their local 
environment) are approaching quantum criticality of the 
uniform droplet model. Integrating expression and 
identifying the local droplet energy q = f +r;+Aj — 5x(0), 
we find the relation between droplet energy and coupling 
strength e; ~ exp(27r/ _1 rj/u) where we have defined 

^ 2VA2 + 2VffA 
A + V¥+ v / A 2 + 2 v /2gA' 

We can now switch from integration over local coupling 
constants to integration over energies via the replacement 
/ r P( r i)(' ' ' — * / r P(ti)(' ■ ■ )dei where the distribution 
of droplet energies is a power-law V(e) ~ e a , a' = fa. 
From Eq. (11) we conclude that the "renormalized Grif- 
fiths exponent" a' decreases as the critical point is ap- 
proached. Scaling arguments show that such a distri- 
bution leads to strong Griffiths behavior in the isolated 
droplet case for a' < 1. 

For the case of interacting droplets, however, integra- 
tion of expression (j^J produces x( w ) — x(Q) ~ ~M" _1 + 
0(M) and non-Ohmic dissipation obtains for a' < 2. As 
argued previously this qualitatively changes the dy- 
namics of the droplets which then find themselves above 
their lower critical dimension. Sufficiently large droplets 
(fi > r c ) can therefore order (freeze) in imaginary time 
before the Griffiths phase is reached, and this new phase 
has been dubbed the "cluster glass" . 

IRFP of isolated droplets at quantum criticality- It 
should be pointed out that our IIL expression for 
a' underscores the close analogy between the present 
study and the physics of infinite-randomness fixed points 
(IRFP). Most remarkably, from Eq. (11) we conclude 
that, within such IIL calculation, the "renormalized 
Griffiths exponent" a' vanishes precisely at the uniform 
droplet quantum critical point as 

a 1 ~ A 1 / 2 . (10) 

This result should be contrasted with the theory of nonin- 
teracting droplets @| where the "bare" Griffiths exponent 
a is a non-universal function of parameters. 

The dynamical critical exponent governing an IRFP is 
commonly defined [l9[ as z' = d/a' (d is the spatial di- 
mensionality) , diverging right at the transition. Equation 
provides us with an explicit expression for z' which 



captures both the quantum Griffiths phase and the be- 
havior reminiscent of an IRFP at A = 0. Of course, this 
expression is not fully self-consistent, and will not be ap- 
propriate at the lowest temperatures where the system 
freezes into the cluster glass state. However, quantum 
Griffiths and IRFP behavior obtained from the IIL ap- 
proximation is expected to hold at temperatures above 
the cluster glass transition, as anticipated previously [f| . 

Fluctuation induced first-order transition. — In the 
paramagnetic phase, it is clear that the leading order 
behavior of x( w ) is linear in lo and it is not difficult 
to show that the slope m x = d\/ 9lo\ u= o determines 
the renormalized exponent a' = fa through the factor 
/ = (1 — gm x ) _1 , as can be checked at the IIL level. This 
relationship can be further illuminated by calculating m x 
self-consistently from equation ©: 

m x = -2(1 - £m x )F(0). (11) 

Solving for the slope, we find the relationship 

/ = (1 - -gm^- 1 = (1 - 2gxHd)). (12) 

What is striking about this result is that the term in 
parentheses on the right-hand side of (|12p is the ex- 
act analogue of the spin glass stability criterion Asg for 
droplets with distributed site-energies [18j; in this case, 
the stability criterion is Asg = 1 — 2gx 2 (0). The fully 
self-consistent value of a' is then a' = fa = Asg«, im- 
plying that a' = at the transition determined by the 
stability of the CG phase. On the other hand, approach- 
ing from the PMFL side, the transition is still given by 
a' = 2 where non-Ohmic dissipation obtains. So, unlike 
the uniform droplet case, the stability criteria for the 
PMFL and the CG do not coincide at a unique point, 
implying that the quantum phase transition in the dis- 
tributed droplet case is not a conventional second-order 
transition. 

To gain further insight into the nature of the transition, 
we recast the problem as an eigenvalue analysis of the 
constitutive free energy for which the solution to © is a 
minimum. The relevant contribution to the free energy 
can be expressed near the minimum as 

T x = I (xH-Xo(w))r x (£J,w')(x(w')-Xo(w')) (13) 

so that the eigenvalue determining the stability of the 
solution xo satisfies 

L,M") -XoM)r x (a;X)(x(a/) - Xo(^)) 

L(xH -xoM) 2 

(14) 

To implement this iteratively, we can approximate the 
difference (x — Xo) ~ (x ~ X^): where x is 
the value of x a t the nth iteration step. To deter- 
mine how T x operates on x^? note that the free en- 
ergy functional is constructed such that ST jb\ = 
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FIG. 1: Upper panel: The boundary corresponding to A x = 
is shown for finite temperatures. The line is a fit to the data 
as discussed in the text. Lower panel: The stability criterion 
Asg calculated at T x . The line is a guide to the eye. In both 
panels, open squares correspond to first-order transitions, full 
circles to second-order transitions. The tricritical point is 
indicated by the solid star. 



yields the self-consistency condition (O for x- Thus, 
J duj'T x {uj,Lj')x {n) {uj') = X (n) M - x (n+1) M. We can 
now rewrite the denning relation (fT4")) for the eigenvalue 
valid at each iteration step 



AgO = i _ 



/dc(x ( "- 1) -x (n) )(x (n) -x ( " +1) ) 

fdw(x ( ™ _1) -x (n) Y 



(15) 



and proceed by iterating numerically. The instability of 
the PMFL is then given by A x = at n — » oo. 

The results of this procedure at finite temperature are 
shown in the upper panel of Figure Q] Plotted in both 
full and empty symbols is the temperature T x at which 
A x = for given values of the tuning parameter f. Plot- 
ted in the lower panel is the value of Asg at T x (f). We 
can see that A SG > A x = for -0.65 < f < f c (T = 0) 
which indicates that there is a region of coexistence of 
the CG and PMFL phases where A x and Asg are both 
finite and positive (i.e. both solutions are stable). The 
transition for f > —0.65 is therefore of first-order, so the 
open squares represent the spinodal line; the other spin- 
odal corresponds to Asg = 0, A x > 0. For f < —0.65, the 
two criteria merge and the usual continuous transition is 
restored. The boundary given by A x = in this regime is 
then the true phase boundary which we can fit using the 
form log(T c ) ~ (f — f c ) _1 , as shown by the full line in the 
upper panel, reflecting the fact that the critical tempera- 
ture is proportional to the number of frozen droplets @ . 
Thus, open squares represent a discontinuous transition 
for f > —0.65 while the full circles represent a contin- 
uous transition for f < —0.65; the tricritical point at 
f = —0.65 is shown by the solid star. 

In summary, we examined quantum criticality of the 



cluster glass phase which typically emerges as disorder 
is introduced in itinerant magnets. We have shown how 
coupling between droplets enhances the quantum Grif- 
fiths effects, leading to singularly enhanced dissipation 
in the vicinity of the quantum critical point. As a result, 
any system with widely distributed droplet sizes cannot 
display conventional quantum criticality, but instead fea- 
tures a first-order transition at T = 0. In contrast, if 
the effects of the largest droplets are eliminated, con- 
ventional criticality is restored, as we found in a related 
model with a bounded distribution of droplet sizes. Sim- 
ilar behavior is also found at sufficiently high tempera- 
tures, where thermal fluctuations impose a cutoff on the 
Griffiths phase anomalies. This implies that fluctuations 
due to these Griffiths phase anomalies are responsible for 
driving the transition first-order. Our results could be of 
direct relevance to experiments on disordered supercon- 
ductors where singular corrections due to rare regions 
have recently been demonstrated (20j . 
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